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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  to  utilize  the  recently  established  ordering,  receiving,  and 
processing  capabilities  for  TerraS AR-X  along-track  interferometric  synthetic  aperture  radar  (along- 
track  InSAR,  ATI)  data  at  the  University  of  Miami's  Center  for  Southeastern  Tropical  Advanced 
Remote  Sensing  (CSTARS)  for  the  development  of  advanced  ATI  data  products,  such  as  surface 
current  fields  in  narrow  rivers,  and  innovative  higher-order  data  products  characterizing,  for  example, 
the  bathymetry  in  coastal  areas,  properties  of  oceanic  internal  waves,  and  the  ocean  surface  wave 
spectrum,  as  derived  from  amplitude  signatures,  interferometric  phase  signatures,  and  polarimetric 
signatures. 

OBJECTIVES 

Within  the  project  period  of  36  months  (which  started  in  spring  2011),  we  intend  to 

►  acquire  and  analyze  a  comprehensive  set  of  ATI  images  of  rivers,  coastal  areas,  and  open-ocean 
features  such  as  internal  waves  and  tropical  storms; 

►  explore  possibilities  of  obtaining  surface  current  estimates  for  challenging  test  sites,  such  as 
very  narrow  rivers  (width  <100  m),  where  most  of  the  backscattered  power  from  a  fast -moving 
water  surface  may  be  mapped  into  pixels  on  land; 

►  test  techniques  for  two-dimensional  vector  current  field  retrievals  from  combined  TerraS  AR-X 
/  TanDEM-X  along-track  InSAR  data  or  combined  single-instrument  along-track  InSAR  data 
from  ascending  and  descending  overpasses  of  a  test  site; 

►  develop  tools  for  the  generation  of  higher-level  data  products  from  along-track  InSAR-derived 
surface  current  fields,  such  as  bathymetric  maps  of  coastal  waters  and  volume  flow  /  mean 
runoff  estimates  for  rivers  of  uncertain  depth;  and 

►  evaluate  the  potential  of  innovative  InSAR-  and  polarimetry-based  techniques  for  fully  two- 
dimensional  surface  wave  retrievals. 
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These  activities  will  make  significant  contributions  to  an  improved  understanding  and  an  optimal 
exploitation  of  the  potential  of  existing  and  future  spaceborne  SAR  systems  for  U.S.  Navy  operations 
and  for  a  variety  of  research  and  engineering  applications. 

APPROACH 

The  project  consists  of  six  main  Work  Packages,  which  are  defined  as  follows. 

Work  Package  1 :  CSTARS  Upgrades.  When  the  proposal  was  written,  we  were  planning  to  implement 
receiving  and  processing  capabilities  for  TanDEM-X  and  COSMO-4  at  CSTARS  (in  addition  to 
existing  capabilities  for  other  satellites)  and  to  evaluate  the  ATI  capabilities  of  COSMO-4  and 
RADARSAT-2. 

Work  Package  2:  Data  Acquisition  and  Raw  Data  Processing.  This  Work  Package  comprises  the 
collection  and  processing  of  a  variety  of  innovative  data  sets  through  CSTARS,  including  new  ATI 
images  of  several  rivers  of  different  widths,  of  coastal  areas  with  a  strong  bathymetry-modulated  tidal 
flow,  test  areas  with  internal  waves  and  eddies;  dual-beam  ATI  data  from  TerraSAR-X  and  TanDEM- 
X  together  for  vector  current  retrievals;  and  combined  along-track  /  cross-track  InSAR  and 
polarimetric  data  of  rivers  and  ocean  scenes  for  a  variety  of  tests.  Some  data  products  were  already  in 
the  TerraSAR-X  /  TanDEM-X  archives  of  DLR,  acquired  during  the  TerraSAR-X  Dual  Receive 
Antenne  (DRA)  mode  campaign  and  the  TanDEM-X  commissioning  phase  in  2010. 

Work  Package  3:  Advanced  Current  Retrieval  Techniques.  In  this  work  package  we  will  improve  our 
current  retrieval  capabilities  by  testing  current  retrievals  over  new  kinds  of  features,  such  as  oceanic 
internal  waves,  and  by  developing  techniques  for  narrow  rivers  and  for  two-dimensional  vector  current 
retrievals.  Due  to  the  azimuthal  displacement  of  moving  targets  in  SAR  images,  most  of  the  signal 
contributions  from  narrow  water  surfaces  may  be  mapped  into  pixels  on  land.  They  must  be  separated 
from  the  contributions  of  stationary  targets  and  moved  back  to  the  correct  location  in  the  image  to 
pennit  meaningful  current  estimates.  We  will  try  an  approach  that  combines  advanced  interferometric 
processing  of  the  two  original  complex  SAR  images  with  physics-based  statistical  modeling  of  the 
contributions  from  water  and  land  and  our  existing  current  retrieval  procedure  for  the  correction  of 
SAR  imaging  artifacts  and  wave  contributions  over  open  waters.  Two-dimensional  vector  current 
retrievals  may  be  obtained  by  combining  data  from  ascending  and  descending  overpasses,  but  a  more 
promising  development  is  be  based  on  the  idea  of  simultaneous  ATI  data  acquisitions  by  TerraSAR-X 
and  TanDEM-X  with  slightly  different  look  directions.  To  acquire  such  data,  one  of  the  satellites  must 
be  rotated  by  some  angle  for  a  limited  time.  A  demonstration  of  the  dual-beam  ATI  technique  with  an 
airborne  system  was  given  Toporkov  et  al.  [1].  Other  ideas  to  be  attacked  within  this  Work  Package 
are  the  combined  use  of  TerraSAR-X  and  TanDEM-X  to  cover  a  larger  effective  swath  or  to  improve 
the  data  quality  within  the  normal  swath;  multi-baseline  along-track  interferometry;  and  a  comparison 
of  current  fields  retrieved  from  ATI  data  with  results  of  the  Doppler  centroid  method  based  on 
conventional  single-channel  SAR  data  (see  [2]  for  an  example). 

Work  Package  4:  Advanced  Wave  Retrieval  Techniques.  Wave  retrieval  algorithms  for  conventional, 
interferometric,  and  polarimetric  SAR  data  will  be  implemented  and  applied  to  images  acquired  in 
Work  Package  1.  Results  will  be  compared  with  each  other,  with  buoy  data,  and  with  wave  spectra 
from  numerical  wave  model  (WAM)  runs.  Techniques  for  combining  results  of  different  wave  retrieval 
methods  (e.g.  based  on  intensities  and  phases  of  interferometric  data  or  on  intensities  and  polarimetric 
coherences  of  polarimetric  data)  will  be  tested,  and  the  feasibility  of  phase-preserving  wave  field 
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retrievals  and  the  detection  of  freak  waves  will  be  evaluated.  Combining  TerraS  AR-X  and  TanDEM-X 
data  with  a  short  time  lag,  a  retrieval  of  frequency-wavenumber  spectra  may  be  possible.  Conclusions 
regarding  optimal  wave  retrieval  strategies  and  desirable  instrument  and  algorithm  improvements  for 
further  improved  wave  retrieval  capabilities  will  be  drawn. 

Work  Package  5:  Higher-Level  Data  Products.  The  retrieval  of  bathymetric  maps  from  ATI-derived 
current  fields,  requiring  only  a  small  number  of  known  depths  for  deriving  a  characteristic  current- 
depth  relationship  for  each  test  area,  has  been  demonstrated  in  [3].  This  existing  algorithm  will  be 
developed  a  little  further  in  such  a  way  that  it  can  be  applied  to  new  test  sites  without  case-specific 
optimization  and  that  phases  and  intensities  of  ATI  images  can  be  analyzed  together.  Furthermore, 
visible  wave  refraction  patterns  will  be  included  in  the  analysis.  Here  we  need  to  account  for  the 
dependence  of  wave  signatures  in  SAR  images  on  the  imaging  geometry  and  for  effects  of  SAR 
imaging  artifacts  such  as  velocity  bunching.  In  addition  to  wave  refraction  due  to  changes  in  the  water 
depth  and  the  corresponding  wave  velocity  according  to  the  shallow-water  dispersion  relation,  some 
contribution  to  the  modulation  of  wavelengths  and  directions  of  propagation  may  result  from  wave- 
current  interaction.  We  think  this  contribution  can  be  identified  and  taken  into  account  in  the  data 
interpretation  if  wave  trains  propagating  in  different  directions  are  considered.  Furthermore,  wave 
patterns  in  SAR  images  acquired  at  different  tidal  phases  will  be  subject  to  different  tidal  current 
fields,  but  a  constant  bathymetry.  At  the  end  of  the  data  processing  chain,  consistent  bathymetric  maps 
should  be  obtained  from  all  images.  In  addition  to  the  bathymetry  work,  an  algorithm  for  river  volume 
flow  estimates  will  be  implemented,  based  on  methods  discussed  in  [4]  and  [5],  and  techniques  for 
separating  tidal  contributions  from  the  mean  river  runoff  will  be  tested  on  the  basis  of  available  long 
time  series  of  ATI  images  of  rivers.  Furthennore,  the  use  of  along-track  and  cross-track  InSAR  data 
from  TerraS  AR-X  and  TanDEM-X  for  combined  surface  curren  and  water  level  measurements  and 
corresponding  improved  volume  flow  estimates  will  be  tested. 

Work  Package  6:  Dissemination.  Results  of  the  project  will  be  disseminated  through  a  project  website, 
presentations  at  international  conferences,  and  publications  in  peer-reviewed  journals.  To  ensure  an 
efficient  collaboration  and  data  exchange  with  existing  ONR  program  teams  and  other  potential 
partners  and  users,  we  will  contact  these  colleagues  in  an  early  stage  of  the  project,  inform  them  about 
our  objectives,  plans,  and  needs,  and  establish  sustainable  communication  links. 

WORK  COMPLETED 

Work  Package  1:  CSTARS  Upgrades.  The  implementation  of  new  receiving  capabilities  at  CSTARS 
was  almost  completed  before  this  project  started.  Software  upgrades  for  a  complete  processing  of 
TerraSAR-X  and  TanDEM-X  ATI  products  at  CSTARS  are  expected  to  be  completed  until  the  end  of 
2011.  The  final  version  of  COSMO-4  does  not  have  ATI  capabilities,  thus  there  is  nothing  to  be  tested 
with  COSMO  satellites  within  this  project.  Regarding  the  availability  of  RADARSAT-2  ATI  data  for 
testing  we  are  more  skeptical  than  we  were  when  the  proposal  was  written,  but  we  know  that  the 
technical  capability  to  run  RADARSAT-2  in  a  split-antenna  mode  like  TerraSAR-X  exists. 

Work  Package  2:  Data  Acquisition  and  Raw  Data  Processing.  In  201 1  we  have  not  acquired  any  new 
ATI  data  so  far,  but  we  have  received  a  couple  of  interesting  datasets  from  the  TerraSAR-X  DRA 
mode  campaign  in  2010.  This  includes  two  ATI  images  of  internal  waves  at  Dongsha  (South  China 
Sea).  We  also  have  TerraSAR-X  ATI  images  that  show  clear  phase  signatures  of  surface  wave 
motions,  and  we  are  currently  negotiating  to  obtain  TerraSAR-X  /  TanDEM-X  image  pairs  from  the 
TanDEM-X  commissioning  phase.  We  will  begin  to  acquire  new  TerraSAR-X  ATI  data  in  Aperture 
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Switching  (AS)  mode  in  fall  /  winter  2011  /  2012,  and  to  downlink  TerraSAR-X  ATI  datasets  to 
CSTARS  for  the  first  time. 

Work  Package  3:  Advanced  Current  Retrieval  Techniques.  Aside  from  some  software  improvements, 
the  main  achievement  in  this  task  so  far  is  the  successful  retrieval  of  currents  over  an  internal  wave  at 
Dongsha  from  a  TerraSAR-X  ATI  image.  To  our  knowledge,  this  has  been  done  with  satellite  data  for 
the  first  time.  A  manuscript  for  a  peer-reviewed  journal  is  in  preparation.  Some  of  the  existing  results 
are  shown  in  the  figures  at  the  end  of  this  report. 

Work  Package  4:  Advanced  Wave  Retrieval  Techniques.  We  do  have  polarimatric  images  of  waves 
from  the  TerraSAR-X  DRA  mode  campaign,  but  we  have  not  started  to  work  on  this  work  package 
within  the  first  months  of  the  project. 

Work  Package  5:  Higher-Level  Data  Products.  We  have  started  to  analyze  wave  patterns  in  a 
TerraSAR-X  spotlight  SAR  image  to  test  our  ideas  about  bathymetry  retrievals  from  wave  refraction 
patterns.  First  results  are  promising,  but  the  work  package  is  still  in  a  very  early  stage.  The  analysis  of 
ATI-derived  current  fields  for  bathymetry  retrievals  has  not  even  started.  Suitable  images  will  be 
acquired  in  2012. 

Work  Package  6:  Dissemination.  We  are  planning  to  submit  a  manuscript  about  the  first  internal  wave 
results  within  the  next  few  months;  possibly  to  a  high-impact  journal  such  as  Science  or  Nature.  We 
will  also  present  results  of  this  projects  at  international  conferences  in  2012,  such  as  AGU  Ocean 
Sciences,  ESA  SeaSAR,  and  IGARSS. 

RESULTS 

The  results  obtained  so  far  indicate  that  it  is  feasible  to  obtain  direct  measurements  of  the  current 
variations  over  an  oceanic  internal  wave  from  TerraSAR-X  ATI  data.  This  pennits  much  more  robust 
estimates  of  internal  wave  amplitudes  and,  thus,  pycnocline  depths  etc.  than  methods  based  on  radar 
intensity  signatures,  since  ATI  signatures  are  much  less  difficult  to  model  and  less  sensitive  to 
uncertainties  in  wind  estimates  (discussed  for  intensity  signatures  in  [6]).  We  are  currently  trying  to 
obtain  in-situ  data  for  validation  of  our  results  and  to  acquire  additional  ATI  images  of  internal  waves. 

IMPACT/APPLICATIONS 

Our  new  algorithms  will  permit  current  and  wave  retrievals  under  conditions  where  they  have  not  been 
possible  before  (e.g.  current  retrievals  in  narrow  rivers,  wave  retrievals  in  regions  with  dominant 
azimuth-traveling  waves)  and  /  or  with  significantly  improved  accuracy.  The  algorithms  for  higher- 
order  products  will  facilitate  the  use  of  SAR  data  for  a  variety  of  applications  and  lead  to  a  better 
utilization  of  existing  spacebome  SAR  capabilities  for  routine  applications. 

RELATED  PROJECTS 

All  ONR-funded  research  projects  dealing  with  remote  sensing  of  rivers,  tidal  flats,  and  current 
features  and  wave  spectra  in  the  open  ocean,  as  well  as  related  U.S.  Navy  op-erations,  can  benefit  from 
this  project.  The  project  has  particularly  strong  connections  to  the  ongoing  "Tidal  Flats  Dynamics" 
DRI,  in  which  CSTARS  and  the  PI  of  this  project  are  directly  involved. 
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Fig.  1:  TerraSAR-X  ATI  image  of  internal  waves  at  Dongsha  (South  China  Sea),  averaged  on  a 
100  m  x  100  m  grid,  trimmed  to  approx.  29  km  x  120  km.  The  left-hand  part  of  the  figure  shows 
SAR  image  intensities  with  a  logarithmic  grayscale.  The  right-hand  part  shows  the  same  with 
superimposed  ATI-derived  line-of-sight  Doppler  velocity  variations  in  color.  Flight  direction  is  191° 
(here:  top  to  bottom),  look  direction  281°  (here:  right  to  left).  Color  scale:  Red  =  0.5  m/s  to  left,  Blue 
=  0.5  m/s  to  right.  White  frames  indicate  areas  for  which  intensity  and  velocity  profiles  have  been 
derived.  Note  that  large  image  intensities  over  the  internal  waves  are  strongly  correlated  with  large 

Doppler  velocities  in  the  propagation  direction. 
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Fig.  2:  Preliminary  result  of  simulations  of  the  internal  wave  signature  at  the  bottom  of  the  image. 
The  diagram  on  the  left  shows  the  surface  current  field  of  a  possible  internal  wave  in  the  test  area, 
parameterized  such  that  good  agreement  between  ATI-derived  and  numerically  simulated  Doppler 
velocities  is  obtained.  The  center  diagram  shows  the  model  result  as  a  bold  line  and  the  measured 
Doppler  velocity  profile  as  derived  from  Fig.  1  as  a  thin  line.  The  agreement  is  good  for  a  peak 
surface  current  on  the  order  of  0.9  m/s  (approx.  0. 7  m/s  in  the  Doppler  velocity).  The  diagram  on  the 
right  shows  that  the  strong  image  intensity  modulation  of  this  signature  cannot  be  reproduced  by 
our  numerical  SAR /ATI  imaging  model;  the  same  has  been  found  in  other  applications  of  the 
model  to  scenarios  with  very  strong  image  intensity  signatures. 


Fig.  3:  Same  as  Fig.  2  for  the  second  white  box  from  the  top  in  Fig.  1.  In  this  case  we  obtain  good 
agreement  of  ATI-derived  and  simulated  Doppler  and  image  intensity  signatures.  Here  we  have 
simulated  the  first  two  solitons  of  the  internal  wave  train  to  obtain  a  realistic  current  gradient 

behind  the  leading  soliton. 
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